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Abstract: Railways built in cold, snowy, and lightly populated areas are subjected to wind and snow 
disasters. In this study, we selected a snow hazard prevention and control section of the Altay-Zhundong 
Railway in Xinjiang Uygur Autonomous Region of China as the research object. We investigated the 
deposited snowfall variation characteristics on the two sides and in the embankment pavement area of 
snow fences with different porosities, fence heights, and arrangement distances using single-factor tests 
and orthogonal tests based on global atmospheric reanalysis climate data, field survey data, and a 
multi-phase flow analysis model. The results showed significant differences in the characteristics of snow 
cover distribution and snow cover thickness between the embankment and the cutting in the absence of 
snow protection measures. The maximum snow cover thickness of the embankment pavement decreased 
by 12.6% relative to the cutting pavement. The snow cover thickness of the embankment exhibited an 
increasing trend from windward shoulder to leeward shoulder, whereas the snow cover thickness of the 
cutting presented a declining trend from windward shoulder to leeward toe. In the collaborative 
prevention and control of snow fences and embankments, the three factors can be ranked in terms of 
their sensitivity to deposited snowfall within the influence scope of snow fences as follows: fence 
height>arrangement distance>porosity. At the same time, fence height yielded a significant relationship 
for the influence scope of snow fences (P<0.05). The three factors can also be ranked in terms of their 
sensitivity to deposited snowfall on the pavement as follows: porosity>fence height>arrangement 
distance. For the embankment protection of the Altay-Zhundong Railway against wind and snow, snow 
fence with a porosity of 75%, a fence height of 4.8 m, and an arrangement distance from the 
embankment of 60 m produced the best snow control effect. By revealing the characteristics of snow 
cover distribution along railway subgrade and the protective effect of snow fences, this study provides 
valuable references for the engineering applications of railway construction in areas prone to wind and 
snow disasters. 
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1 Introduction 


Snowdrifts experience creep, saltation, and suspension in the flow field and redistribute snowfall 
regularly, causing snow depths in some areas to reach 3—8 times the average snow depth (Xia et 
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al., 2017). However, the established railway lines of traffic network construction projects are 
rarely changed simply due to the harsh natural environment, making them vulnerable to wind and 
snow disasters. During snow melting, the mixture of ice, snow, and water infiltrates the railway 
subgrade, resulting in settlement (or collapse) and traffic interruption (Lu et al., 2013). Owing to 
the complexity of railway track structures, light snow cover often limits the execution of train 
traction, braking, and bogie systems, while heavy snow cover may even lead to out-of-service 
maintenance (Mott et al., 2018; Wang et al., 2018). In order to maintain the railways' working 
efficiency and stable operation and reduce late-stage input costs, it is necessary to identify the 
characteristics of snow cover distribution along railways and feasible to take appropriate 
protective measures around railways. 

Field surveys, wind tunnel tests, and computational fluid dynamics (CFD) numerical 
simulation are different methods to study the characteristics of snow cover distribution and 
protective measures. Field research provides researchers with information about site geological 
conditions, enabling them to evaluate the actual effects of protective measures (Beria et al., 2018; 
Liu et al., 2021). In this regard, field surveys, as a direct and effective method, obtain wind speed, 
snow thickness, start-up conditions, and other information in the field for laboratory simulation 
(Gordon and Taylor, 2008; Loeb and Kennedy, 2021). However, it is challenging to conduct 
precise quantitative research for the railway lines passing through areas prone to wind and snow 
disasters or sparsely populated Gobi areas, mainly due to the large accidental errors and low 
reproducibility of measurement and the presence of many uncontrollable external factors 
(Tominaga, 2018). Researchers have widely recognized wind tunnel tests due to the advantages of 
controllable environmental variables and systematic parametric measurement (Cao et al., 2019; 
Kumar et al., 2020). Wind tunnel tests allow scaling building entities by adopting granular 
materials with different properties and performing calculations according to the similarity 
criterion (Zhou et al., 2019). However, wind tunnel tests are inadequate, expensive, and adaptable 
to the similarity criterion to varying degrees (Tominaga et al., 2011; Gao et al., 2022). 

In recent years, with the optimization of simulation calculation and theoretical models, many 
researchers have widely adopted CFD numerical simulation to study snowdrifts (Zhou et al., 2021). 
Different forms of snow fences can be erected to effectively alleviate the impact of wind and snow 
disasters on railway operations (McClure et al., 2017). The control effects of snow fences are related 
to structural factors, such as porosity, arrangement distance, fence height, and arrangement form of 
multiple rows of snow fences (Kumar, 2015), and environmental factors like wind direction, wind 
frequency, snow thickness, topography, geomorphology, elevation, and complex wind and snow 
conditions. Thus, it is difficult to establish specific design standards or take all-purpose protective 
measures (Potac et al., 2018). Basnet et al. (2014) highlighted that the snow accumulation in the 
bottom gap of snow fence is often not significant, and this design would be effective for areas with 
severe snow-drifting hazards. Liu et al. (2016) examined the distribution of snow deposition and 
erosion on both sides of temporary snow fences considering various porosities, fence heights, and 
bottom holes. This study further established a simplified relationship between fence characteristics 
and the site of snow deposition. Prokop and Procter (2016) modified the existing snow fence 
structure and applied it in the field based on the topography and wind data. Currently, there is a lack 
of theoretical support for predicting the characteristics of snow cover distribution along railway 
subgrade, and the protective effects of different forms of snow fences remain to be clarified. In 
order to improve prediction precision, it is imperative to explore the inducements of snow cover 
distribution and study the collaborative prevention and control effect of snow fences and subgrade. 

Some sections along the Altay—Zhundong Railway are plagued by snow damage, which 
seriously threatens the railroad's safe operation. The impact of snow damage on the railway line 
must be observed and monitored. Snow damage prevention and control is the key to the railroad's 
safe operation. Due to the uninhabited railroad line during the snowfall period and the lack of 
meteorological information, it is difficult to conduct a study on the snow damage prevention and 
control for the Altay—Zhundong Railway. Based on the erosion—accumulation model (Wang et al., 
2019) and the multi-phase turbulent flow theory (Cheng et al., 2016; Zhang et al., 2021), in this 
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study, we used global atmospheric reanalysis (ERAS, where ERA5 is the fifth generation 
ECMWFE (The European Center for Medium-Range Weather Forecasts) atmospheric reanalysis of 
the global climate covering the period from January 1950 to present) climate data to validate and 
parameterize the wind speed, wind frequency, snow density, snow thickness, and other data 
acquired from field surveys, thereby supporting laboratory simulation tests. On this basis, we 
analyzed the collaborative prevention and control effect of snow fences and subgrade along the 
Altay—Zhundong Railway from the perspectives of porosity, fence height, and arrangement 
distance and designed orthogonal tests to explore the sensitivity of each factor to different areas 
and related prevention and control effects. This study provides theoretical guidance and technical 
support for predicting snow cover distribution and constructing snow protection facilities along 
railways in cold, windy, and snowy areas. 


2 Materials and methods 


2.1 Study area 


2.1.1 Regional geographical characteristics 

The Altay—Zhundong Railway (88°05'24"—89°33'48"E, 44°44'24"47°43'12"N), also known as 
Altay—Fuyun—Zhundong Railway (Fig. 1), is located in the northern Xinjiang Uygur Autonomous 
Region, China. It starts from the Altay Railway Station in the north, ends at the Zhundong North 
Railway Station in the south, and connects with Urumqi—Jiangjunmiao Railway at the north end. 
The Altay—Zhundong Railway, with a total length of 420.4 km, is a key railway of the Northern 
Xinjiang Railway Loop. This study mainly investigates the wind and snow disasters threatening 
the Altay—Fuyun section. 


Fig. 1 Overview of the Altay—Zhundong Railway 


Built at the southern edge of the Altay Mountains and north of the Junggar Basin, the 
Altay—Zhundong Railway belongs to the hinterland of Eurasia and has a typical cold temperate 
continental climate. The Altay-Zhundong Railway is persistently subjected to snow hazards, 
severely impacting local traffic operation due to the cold winter, frequent rains, thick snow cover, 
and long snow cover duration (Fig. 2). This makes it necessary to erect snow protection facilities 
along the railway to prevent it from being covered by snow. However, existing snow protection 
facilities have uneven effects, and the railway still suffers from damage by snow hazards. 
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cm 


Fig. 2 Snow hazards along the Altay—-Zhundong Railway. (a), snow cover along the railway line; (b), snow 
cover alongside slopes; (c), snow cover along snow fences. 


2.1.2 Regional snow source conditions 

Abundant snow sources are necessary for the formation of snow disasters. In order to identify 
regional snow cover patterns, in this study, we selected a site (88°41'24"E, 47°21'00"N) along the 
railway as the research focus. Global atmospheric reanalysis (ERA5) climate data were taken as 
the data source to acquire the study area's snow depth and snow density data in 2018, 2019, and 
2020. According to a statistical analysis of data, the annual average snow density of the study area 
in 2018, 2019, and 2020 was 132.26, 134.94, and 130.99 kg/m3, respectively, and the average 
snow density of the study area in the three years was 132.73 kg/m>. The main forms of snowfall 
in the study area were fresh snow (wet snow) and fine snow (generally less than 0.50 mm in grain 
size). The threshold wind speed was 3.7—4.3 m/s, with a median of 4.0 m/s. The threshold wind 
speed increased linearly with the square root of snow grain size at a temperature lower than —6°C, 
as expressed in the following equation (Ma and Huang, 2006). 


V,=344+1.5VD, (1) 


where JV; is the threshold wind speed (m/s); and D is snow grain size (mm). 

The range of snow grain size was calculated as 0.04—0.36 mm, with a median of 0.20 mm. The 
number of snow-covered days in the study area in 2018, 2019, and 2020 was 162, 161, and 170 d, 
respectively, lasting from October to April of the following year. The number of days with a snow 
thickness above 10 cm in 2018, 2019, and 2020 was 50, 111, and 56 d, respectively. The 
maximum snow thickness of the study area in the three years was 44 cm (Fig. 3). Generally, the 
study area is covered by snow for about half of the year, and snow thickness exceeds 10 cm in 
about 2-3 months of a year. The thick snow cover and long snow cover duration provide 
favorable conditions for forming snow disasters. 
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Fig. 3 Time-history diagram of snow thickness in the study area from 2018 to 2020 
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2.1.3 Regional wind conditions 


In addition to snow sources and geomorphological conditions, strong winds also contribute to the 
formation of wind and snow disasters. The wind data of the study area in the three years were 
derived from ERAS. The statistical analysis of data revealed that the west-northwest (WNW) 
azimuth frequencies of the study area in 2018, 2019, and 2020 reached 53.45%, 54.55%, and 
57.29%, respectively, and that the east-northeast (ENE) azimuth frequencies in the winter of 
2018, 2019, and 2020 reached 73.49%, 89.55%, and 81.51%, respectively. According to wind 
speed statistics, the time when the wind speed exceeded the threshold wind speed in the study 
area in 2018, 2019, and 2020 accounted for 16.56% of the total time of the whole year. The 
dominant direction was WNW; the average and maximum wind speeds were 5.56 and 12.57 m/s, 
respectively. The time when the wind speed exceeded the threshold wind speed in the winter 
accounted for 9.55% of the total time of the winter. The dominant direction was ENE. The 
average and maximum wind speeds were 4.32 and 9.38 m/s, respectively. Overall, the maximum 
WNW azimuth frequency of the study area in the three years was 57.29%, and WNW was the 
prevailing wind direction in the whole year. The maximum ENE azimuth frequency in the winter 
was 89.55%, and ENE was the prevailing wind direction. The prevailing wind direction in the 
winter differed greatly from that in the year (Fig. 4). The study area showed the characteristics of 
intermittent snow transport and stable snow transport direction. The occurrence of a wind and 
snow disaster necessarily means that there has been snowfall, so winter data are more precise and 
reasonable. It explains the justification for selecting ENE as the prevailing wind direction in this 
study. In this study, ERA5 data were used to obtain the wind direction, wind speed, and real-time 
snow thickness in the study area. Considering the snow thickness and wind speed simultaneously, 
the wind speed was determined when the snow thickness was large (the necessary conditions for 
forming a wind-blowing snow disaster). Combined with the actual layout of the project, the snow 
thickness obtained from ERA5 data helped to determine wind speed. The obtained wind speed 
was further applied to the inlet wind speed during the numerical simulation. 
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Fig. 4 Rose plot of the frequency of wind directions. (a), 2018; (b), 2019; (c), 2020; (d), in the winter of 2018; 
(e), in the winter of 2019; (f), in the winter of 2020. MWS, max wind speed; AWS, average wind speed; N, north; 
NNE, north-northeast; NE, northeast; ENE, east-northeast; E, east; ESE, east-southeast; SE, southeast; SSE, 
south-southeast; S, south; SSW, south-southwest; SW, southwest; WSW, west-southwest; W, west; WNW, 
west-northwest; NW, northwest; NNW, north-northwest. 
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2.2 Numerical simulation 


2.2.1 Geometric model and grid division 
The buildings on site were reproduced by a scale of 1:1 in Unigraphics NX modeling to ensure 
authenticity and effectiveness. Snow fence model involved multiple parameters, including fence 
height, porosity, and arrangement distance, detailed in part about the design of single-factor and 
orthogonal tests. Modeling was illustrated using a snow fence with a fence height of 4.8 m, a 
porosity of 75%, and an arrangement distance of 20 m. Fence width, embankment height, slope, 
and pavement width were 3.2 m, 2.0 m, 1.0:1.5, and 8.0 m, respectively (Fig. 5). The calculation 
domain size and pavement width of the cutting model and embankment model are kept the same. 
The depth of the cutting is 2.0 m, the embankment height is 2.0 m, and the slope is 1.0:1.5. The 
spatial resolution of ERAS data is about 0.25°. The coordinates along the Altay—Zhundong 
Railway selected in this study are 88°41'24"E and 47°21'00"N. According to the point coordinates 
of the study area, the ERAS data containing the minimum range of known coordinates were 
downloaded, and the other data, except the known coordinates, were discarded. Considering the 
influence of calculation accuracy and efficiency, the computational domain was 150.0 mx32.0 mx 
30.0 m (lengthxwidthxheight). The model is illustrated in Figure 6. 


Pavement width 


Arrangement distance 


Fig. 5 Schematic diagram of the building entity (a) and the computational domain (b) 
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Fig. 6 Schematic diagram of mesh generation for the computational domain 


Unstructured tetrahedral grid division was adopted, and snow fences, embankments, and 
ground grids were densified. Grid independence was validated before calculation to ensure the 
reliability of calculation results. Related grid models were built by adjusting the size parameters 
of snow fences, embankments, and ground grids. The minimum grid size for the model was set as 
0.18, 0.16, 0.15, 0.14, 0.13, 0.12, 0.11, 0.10, 0.09, and 0.08 m. An incoming wind with a constant 
speed of 5.0 m/s was applied in the vertical direction of snow fence. The average wind speed of 
the monitoring surface 10 m behind the entrance of the computational domain was taken as the 
index. The results showed that, with the decrease in grid size, the number of grids multiplied, and 
the wind speed variation rate of the monitoring surface first decreased and then tended to be 
stable (Fig. 7). At the same time, as the number of grids multiplied, calculation time multiplied as 
well. Considering the effects of calculation precision and efficiency, the minimum grid size was 
set as 0.09 m, the number of grids was 3.70x10°, the relative error was about 0.3%, and the 
average grid quality was rated as 0.83, which was excellent. In the simulation, the governing 
equation's residual declined to 107, and the stability of the monitored physical quantity was taken 
as the convergence criteria. It was solved using double-precision solving. 
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Fig.7 Variation rates of wind speed of the monitoring surface under different numbers of grids 


2.2.2 Parameters and boundary conditions 


A multi-phase flow analysis model was built based on CFD numerical simulation. Userdefined 
function (UDF) was used to load the model and optimize boundary conditions. The Euler 
two-fluid method was adopted to calculate the steady state first and the transient state 
afterward. The standard K-epsilon model was adopted. Air was selected as the major phase, and 
snow grains were taken as the minor phase. The default air parameters of the software were 
used as major phase parameters. The average snow density (132.73 kg/m°) was adopted as the 
density of snow grains (minor phase) in the three-year period, and a 0.20 mm grain size was 
adopted as the median. Under initial conditions, the volume fraction of snow grains was set to 
0.05. The variation trend of the entrance wind speed along height was expressed by a 
logarithmic function. The equation of the logarithmic function of the inlet wind speed is as 
follows: 


v 
vo)= hn, (2) 
k Vo 
where v is friction wind speed (m/s); y is the height (m); k is the von Karman coefficient (taken as 
0.40); and yo is the ground roughness (taken as 0.0005 m) (Zhang, 1985). 
The remaining boundary conditions are shown in Table 1. 


Table 1 Setting of the remaining boundary conditions 


Boundary Type 
Outlet Outflow 
Outer side of the computational domain Symmetry 
Ground, subgrade, and snow fence Wall 


2.2.3 Mathematical models and equations 


During snow transport, the effect of snow grain sedimentation should be considered in the 
convective term for the suspension layer. The control equation is given as follows (Kang et al., 


2018): 
ag aU, awh d ə% 
ot ax, Ox, Ox, ‘dx, 


J 


(3) 
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pelts (4) 


where ¢ is the snow concentration; f is time (s); wy is the snowfall speed (m/s); D; is the turbulent 
diffusion coefficient; U; is the velocity vector component of wind; x; is the vector component of 
spatial coordinates; x3 is the vertical component of spatial coordinate; w is the turbulent kinematic 
viscosity; and Sc; is the turbulent Schmidt number. 

The model proposed by Naaim et al. (1998) was used to evaluate the erosion and deposition of 
snow grains in the simulation, as given in Equations 5 and 6. 


2 2 
qero = Aero (us ~ Ux; ), Us 2 Ust, (5) 
2 2 
Ux, — Ux 
dep = PsW yp Ue S Usg, (6) 
Ux, 


where Aero is the erosion constant, set as —7.0x10%; qero is the erosion flux (kg/(m7-s)); qaep is 
deposition flux (kg/(m7-s)); u, is the wall friction velocity (m/s), calculated by us = 0.02205" , 


where vs denotes the reference height wind speed (m/s; Beyers et al., 2004); u»: is the threshold 
friction velocity of snow grains (m/s), set as 0.15 m/s; and @, is the mass concentration of snow 
(kg/m3), expressed as @=p,f (where ps denotes snow grain density (kg/m*) and f denotes 
snow-phase volume fraction). 

For the unification of the calculated quantities in simulation, qtotaı (kg/(m*-s)) was used as a 
substitute for qero and dep, 1-€., Ytotal is the amount of snow cover on the ground that can be 
considered a product of the erosion and deposition motion of snow grains. The amount of snow 
cover on the ground produced by snow grain motion within time Aż (s) is defined as Ah (m), 
which can be expressed as follows: 


Ah=-— ditotal At : (7) 
Ps¥ 


where y is the maximum volume fraction of snow cover, set as 0.62. 
2.3 Test design 


2.3.1 Design of single-factor tests 

In general, fence height should be set at the moderate level of 2.0—6.0 m, as excessively low snow 
fences produce a poor snow resistance effect, while excessively high ones have compromised 
stability and construction costs (Sanudo-Fontaneda et al., 2011). To increase snow cover on both 
sides of snow fence and facilitate the analysis of protective effect, we set the invariant fence 
height, arrangement distance, and porosity as 4.8 m, 80 m, and 75%, respectively. Computational 
variables and working conditions are provided in Table 2. 


Table 2 Computational variables and working conditions 


Working condition Porosity (%) Fence height (m) Arrangement distance (m) 
1 0, 25, 50, 75 4.8 80 
2 75 2.1, 3.0, 3.9, 4.8 80 
3 75 4.8 20, 40, 60, 80 


2.3.2 Design of orthogonal tests 

We designed L16 (4°) orthogonal tests according to three factors (including fence height, 
arrangement distance, and porosity) in single-factor tests. Orthogonal conditions and orthogonal 
results are provided in Table 3. The area within the influence scope of the off-embankment fence 
was designated as Area I. The scope was set as from 20 m in front of snow fence to the 
embankment behind snow fence. The embankment pavement area was designated as Area II. 
Variance and range analysis methods were used to analyze the factors affecting the deposited 
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snowfall on both sides of snow fence and the embankment pavement area. 


Table 3 Orthogonal conditions and orthogonal results 


Average volume fraction of snow grains 


NO. Porosity (%) Fence height (m) Arrangement distance (m) 

Area I Area II 

1 0 2.1 20 0.39273 0.12124 
2 0 3.0 40 0.27658 0.11452 
3 0 3.9 60 0.35911 0.12174 
4 0 4.8 80 0.44215 0.10029 
5 25 2.1 40 0.36019 0.12694 
6 25 3.0 20 0.33829 0.10029 
7 25 39 80 0.39491 0.14488 
8 25 4.8 60 0.48909 0.10821 
9 50 2.1 60 0.29172 0.14166 
10 50 3.0 80 0.30679 0.13402 
11 50 3.9 20 0.50527 0.10015 
12 50 4.8 40 0.49930 0.08282 
13 75 2.1 80 0.29993 0.15682 
14 7 3.0 60 0.39738 0.13629 
15 75 3.9 40 0.40315 0.10301 
16 75 4.8 20 0.50140 0.29079 


Note: Area I and Area II represent the area within the influence scope of the off-embankment fence and the embankment 
pavement area, respectively. 


3 Results 


3.1 Comparative analysis between the embankment and the cutting 


The calculation of the embankment and the cutting was performed in the absence of a snow fence, 
and it was found that there were significant differences in the characteristics and thickness of 
snow cover distribution between the embankment and the cutting. By comparing their maximum 
snow cover thickness points, we found that the maximum dimensionless snow cover thicknesses 
of the embankment and the cutting were 0.153 and 0.175, respectively. The maximum snow cover 
thickness of the embankment pavement decreased by 12.6% relative to the cutting pavement. The 
snow cover thickness of the embankment presented a gradually rising trend from windward 
shoulder to leeward shoulder, whereas the snow cover thickness of the cutting presented a 
gradually declining trend from windward shoulder to leeward toe (Fig. 8). Due to the blocking 
effect of the embankment on wind speed, wind speed formed a vortex at the windward slope toe 
of the embankment. The flow direction of the vortex was negatively correlated with the wind 
speed of the upper layer, which reduced wind speed in the study area and weakened the ability of 
wind to transport snow grains. As a result, snow thickness presented a transient rising trend near 
the windward slope toe of the embankment. In the windward slope toe-shoulder stage of the 
embankment, the snow cover thickness on the ground dropped abruptly because the wind speed of 
the upper layer was not blocked or interfered by the embankment blocking. In contrast, the 
deposited snowfall variation on the slope surface depended on the snow-carrying capacity of 
upper-layer air. In the pavement stage, snow grains on the ground were carried by wind speed. 
Hence, friction velocity exceeded the threshold, and snow cover thickness increased rapidly along 
the wind direction. In conclusion, the embankment itself has a protective effect against wind and 
snow disasters. It is suggested that embankments be preferentially built in areas prone to wind and 
snow disasters. 
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Fig. 8 Comparative analysis of dimensionless snow cover thickness between embankment height (a) and cutting 
depth (b) 


3.2 Flow field distribution around snow fence 


Figure 9 compares the flow fields with and without snow fences in the Altay-Zhundong Railway. 
It can be seen from Figure 9a that the airflow velocity was redistributed after the disturbance of 
the embankment, and the airflow deceleration zone and the turbulence zone were formed on the 
windward slope and the leeward slope, respectively. The airflow on the top of the embankment 
was divided into the lower layer, where the airflow velocity near the surface was greatly reduced, 
forming a weak wind zone, and the upper layer, where the airflow velocity was strengthened to 
form a strong wind zone. Since the attenuation of the airflow velocity is the main cause of snow 
deposition, the windward slope and leeward slope foot were the main snow accumulation areas. 

It can be seen from Figure 9b that the reasonable setting of snow fence on the windward side of 
the embankment can effectively change the flow field distribution around the embankment. When 
the airflow was close to anti-snow fence on the windward side of the embankment, the normal 
airflow was stopped, and the airflow velocity on the windward side of the anti-snow fence was 
decreased. The airflow increased rapidly and reached the top of snow fence as it passed through 
the snow barrier. As a result of the reverse pressure gradient, a reverse flow area formed on the 
leeward side of snow fence, which led to the settlement of snow near the leeward side of snow 
fence. As a consequence, the snow carried by the wind was greatly reduced, effectively reducing 
the settlement of snow on the windward slope of the embankment. For the leeward side of the 
embankment, the range of the leeward slope turbulence zone was greatly reduced, effectively 
reducing the settlement of snow at the foot of the leeward slope. For the embankment pavement, 
the reasonable setting of snow fence increased the wind speed of the pavement, thus increasing 
the ability of the wind to carry snow and effectively reducing the snow settlement on the 
embankment pavement. In summary, setting snow fences to change the flow field distribution 
around the embankment can affect and control the number and settlement trend of snow particles 
transported on the subgrade's windward slope, road surface, and leeward slope. 


3.3 Analysis of single-factor tests 


3.3.1 Porosity 

Figure 10 shows the deposited snowfalls under different porosities. Only black lines and white 
balls are used to display some streamlines and snow grains to reduce overlapping between cloud 
atlases. It can be seen that the snow cover distribution patterns and streamlined vortex forms in 
front of and behind snow fence varied greatly with changing porosity. 

Under porosities of 0% (impervious snow fence), 25%, 50%, and 75%, the average volume 
fractions of snow grains on the two sides of snow fence (within 10 m scope in front of and behind 
snow fence) were determined as 0.523, 0.502, 0.475, and 0.530, respectively, while those in front 
of the embankment (within 10 m scope) were 0.536, 0.532, 0.522, and 0.493, respectively. Also, 
the average volume fractions under porosities of 0%, 25%, 50%, and 75% on the windward side 
slope were 0.191, 0.187, 0.184, and 0.192, respectively, whereas those on the railway pavement 
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Fig. 9 Comparison diagram between the flow field with and without snow fence on the Altay—Zhundong 
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with snow fence, 0 on the x-axis denotes the position of snow fence. The positive and negative values of the 
x-axis represent snow fence's leeward and windward sides, respectively. 
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Fig. 10 Volume fraction of snow grains under different porosities. (a), porosity=0%; (b), porosity=25%; (c), 
porosity=50%; (d), porosity=75%. The value of 0 is the position of snow fence, and the positive and negative 
values represent the leeward and windward sides of snow fence, respectively. 
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were 0.108, 0.111, 0.105, and 0.103, respectively. Under porosities of 0%, 25%, and 50%, the 
vortex formed in front of snow fence was transformed into a jet on snow fence under the action of 
pores, and friction velocity increased suddenly, causing snow to be blown away from both sides 
of snow fence. However, when porosity was 75%, the large porosity produced a weak wind 
resistance capacity, making it difficult to form a high wind speed in pores or a jet. Consequently, 
more snow accumulated on both sides of snow fence. When porosity increased from 0% to 75%, 
the deposited snowfall on both sides of snow fence first decreased and then increased, and the 
inflection point occurred under the porosity of 50%. The deposited snowfall in front of the 
embankment constantly decreased, and the accumulation and transport capacity of the snow cover 
behind snow fence was continuously weakened. The deposited snowfall on the windward side 
slope first decreased and then increased, and the inflection point occurred under the porosity of 
50%. The deposited snowfall on the pavement first increased and then decreased, and the 
inflection point occurred under the porosity of 25%. 

Under the porosity of 75%, snow fence had the strongest snow resistance capacity and the 
largest deposited snowfall on both sides. The snow cover's transport capacity behind snow fence 
was the weakest, and the deposited snowfall in front of the embankment and pavement reached 
their smallest points, respectively. Snow fence produced the best protective effect against 
snowdrifts in this case. 


3.3.2 Fence height 


Figure 11 shows the deposited snowfalls under different fence heights. It can be seen that the 
snow cover distribution patterns and streamlined vortex forms in front of and behind snow fence 
varied greatly with changing fence height. Under 2.1, 3.0, 3.9, and 4.8 m fence heights, the 
average volume fractions of snow grains on the two sides of snow fence were 0.392, 0.443, 0.327, 
and 0.359, respectively, whereas those before the embankment on the pavement were 0.461, 
0.401, 0.626, and 0.601, respectively. The average volume fractions under 2.1, 3.0, 3.9, and 4.8 m 
fence heights were 0.159, 0.150, 0.120, and 0.104, respectively. 

With the change in fence height, the vortex forms before and after snow fence also changed. 
When the fence height was 2.1 m, several small vortexes were formed behind snow fence. These 
small vortexes were far from each other and did not cause a superposition of wind speeds. The 
transport capacity of snow grains was weak, and deposited snowfall concentrated behind snow 
fence. A fence height of 3.0 m led to the formation of mesoscale vortexes behind snow fence and 
the formation small vortexes in front of snow fence. Mesoscale vortexes concentrated in a narrow 
area behind snow fence would not improve the transport capacity of snow grains. As a result, 
snow grains were mainly concentrated behind snow fence, and only a few were deposited in front 
of snow fence. When the fence height was 3.9 m, two relatively large vortexes were formed 
behind snow fence, and a mutually reinforcing trend emerged in the middle area behind snow 
fence, improving the transport capacity of snow grains and expanding the distribution scope of 
snow deposition. When the fence height was 4.8 m, two large vortexes were formed behind snow 
fence, and the transport capacity of snow grains was greatly improved. However, due to the 
reasonable setting of arrangement distance, the improved transport capacity of snow grains would 
not significantly affect the deposited snowfall on the pavement. Considering only the snow cover 
on the pavement, a fence height of 4.8 m produced the largest area of the leeward side snow 
cover, the smallest deposited snowfall on the pavement, and the best protective effect. 


3.3.3 Arrangement distance 


Figure 12 shows the deposited snowfalls under different arrangement distances. It can be seen that 
the snow cover distribution patterns and streamline vortex forms in front of and behind snow 
fence and on pavement varied greatly with changing distance between snow fence and subgrade. 
Under 20, 40, 60, and 80 m arrangement distances, the average volume fractions of snow grains 
on the two sides of snow fence were 0.405, 0.405, 0.386, and 0.344, respectively, those before the 
embankment were 0.619, 0.624, 0.554, and 0.621, respectively, and those on the pavement were 
0.121, 0.099, 0.083, and 0.103, respectively. When the arrangement distance increased from 20 to 
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Fig. 11 Volume fraction of snow grains under different fence heights. (a), fence height=2.1 m; (b), fence 
height=3.0 m; (c), fence height=3.9 m; (d), fence height=4.8 m. The value of 0 is the position of snow fence, and 
the positive and negative values represent the leeward and windward sides of snow fence, respectively. 


80 m, the deposited snowfall on both sides of snow fence continuously decreased with increasing 
distance between snow fence and embankment. The deposited snowfall on pavement continuously 
decreased with increasing arrangement distance from 20 to 60 m, showed an inflection point at 60 
m, and presented a rising trend with increasing arrangement distance from 60 to 80 m. The 
snowfall deposited in front of the embankment was the largest under an arrangement distance of 
40 m, whereas it was the least under an arrangement distance of 60 m. When the arrangement 
distance was 40 m, snow grains began to span across the embankment, and the leeward deposited 
snowfall of the embankment increased significantly. When the arrangement distance decreased 
from 40 to 20 m, the leeward deposited snowfall of the embankment continuously increased. 
When the arrangement distance was 20 m, the vortex scope behind snow fence covered a part of 
the pavement area, and the snow deposition area acted on the pavement. The deposited snowfall 
on the pavement was far greater than that in the absence of a snow fence and that under other 
arrangement distances, in which case snow fence produced no snow protection effect. 

Under a fence height of 4.8 m and a porosity of 75%, the optimal arrangement distance between 
snow fence and embankment was 60 m. In this case, the deposited snowfall in front of the 
embankment and that on the pavement reached their smallest points, respectively, and snow fence 
produced the best protective effect. 
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Fig. 12 Volume fraction of snow grains under different arrangement distances. (a), arrangement distance=20 m; 
(b), arrangement distance=40 m; (c), arrangement distance=60 m; (d), arrangement distance=80 m. The value of 0 
is the position of snow fence, and the positive and negative values represent the leeward and windward sides of 
snow fence, respectively. 


3.4 Analysis of orthogonal tests 


3.4.1 Deposited snowfall within the influence scope of the off-embankment fence 


Table 4 provides the variance analysis of the average volume fraction of snow grains in Area I. 
The difference relationship produced by the three factors (i.e., arrangement distance, fence height, 
and porosity) for the Area I can be solved according to Table 4. The model's coefficient of 
determination (R°) is 0.841, indicating that 84.1% of the variance in Area I was caused by 
arrangement distance, fence height, and porosity. It can be known from the following table that 
fence height produced a significant difference relationship for Area I (P<0.05) and that neither 
arrangement distance nor porosity produced any difference relationship for Area I. 

Table 5 provides the range analysis of the average volume fraction of snow grains in Area I. 
The range values of porosity factor, fence height factor, and arrangement distance factor were 
3.31x107, 1.53x107, and 7.3510, respectively. The three factors can be ranked based on their 
sensitivity to deposited snowfall in Area I as follows: fence height>arrangement 
distance>porosity, i.e., fence height was the main factor affecting the deposited snowfall within 
the influence scope of the off-embankment fence. Porosity and arrangement distance greatly 
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affected the distribution patterns and accumulation locations of snow cover on both sides of snow 
fence but only slightly affected the deposited snowfall in Area I. 


Table 4 Variance analysis of the average volume fraction of snow grains in Area I 


Factor Sum of square Degree of freedom Mean square F P 
Arrangement distance (m) 0.011 3 0.004 1.563 0.293 
Fence height (m) 0.063 3 0.021 8.615 0.014" 
Porosity (%) 0.003 3 0.001 0.409 0.752 

Residual 0.015 6 0.002 - - 


Note: -, no data; *, P<0.05. 


Table 5 Range analysis of the average volume fraction of snow grains in Area I 


Average volume fractions of snow grains Average volume fraction of snow grains 
Level index Porosity Fence height Arrangement Level index Porosity Fence height Arrangement 
factor factor distance factor factor factor distance factor 

i, 0.36764 0.33614 0.43442 k 0.40077 0.41561 0.38433 

ky 0.39562 0.32976 0.38480 ky 0.40046 0.48299 0.36094 


Note: k, k>» k, and k, represent the average value of the sum of the test objectives corresponding to each level, 


respectively, and the values are ranked in descending order under different factors. 


3.4.2 Deposited snowfall at the top surface of the embankment 

Table 6 provides the range analysis of the average volume fraction of snow grains in Area II. The 
range values of porosity factor, fence height factor, and arrangement distance factor were 
5.73x10°, 2.81107, and 2.72%107, respectively. The three factors can be ranked in terms of 
their sensitivity to deposited snowfall in Area II as follows: porosity>fence height>arrangement 
distance. Hence, porosity was the main factor affecting the deposited snowfall in the pavement 
area. Fence height and arrangement distance greatly affected off-embankment snow cover 
distribution and deposited snowfall but only slightly affected the snowfall in Area II. 


Table 6 Range analysis of the average volume fraction of snow grains in Area II 


Average volume fraction of snow grains Average volume fraction of snow grains 
Level - Level 
index Porosity Fence height Arrangement index Porosity Fence height Arrangement 
factor factor distance factor factor factor distance factor 
k 0.11445 0.13667 0.15312 k, 0.11466 0.11744 0.12698 
A 0.12008 0.12128 0.10682 A 0.17173 0.14553 0.13400 


4 Discussion 


Using field surveys to validate remote sensing meteorological data for laboratory simulation tests 
produced a desirable effect in evaluating complex snow disaster prevention and control system 
(Qi, 2018; Gossart et al., 2020; Potsis et al., 2023). ArcGIS, ERAS, and field survey data were 
combined to reveal the study area's geographical location and wind and snow conditions. The data 
obtained from the study area, including snow density, grain size, wind frequency, wind speed, and 
other data, were used for numerical simulation (Tominaga and Stathopoulos, 2020). The 
erosion—accumulation model was used to compare the embankment and the cutting. It was found 
that there were significant differences in the characteristics and thickness of snow cover 
distribution between the embankment and the cutting, and the embankment produced a more 
significant protective effect against snow drifts. Embankments should be preferentially built in 
areas prone to wind and snow disasters without considering the effects of objective factors such as 
economy, space, and laying difficulty. In addition, snow fences help improve the embankment's 
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snow protection effect (Qi et al., 2021). The collaborative prevention and control effect of snow 
fences and subgrade should be considered in the erection of snow fences. Protective measures for 
the embankment and the cutting should be investigated by types with a view to engineering 
practice (Yang et al., 2022). 

This study reveals that, by analyzing the collaborative prevention and control effect of 
pavement and snow fences with different porosities, fence heights, and arrangement distances, 
erecting snow fences significantly improved the snow protection effect of the embankment. A 
parametric analysis of the three factors showed that fence height affects the two sides of snow 
fence more significantly than porosity. The erection of snow fences should consider the spanning 
of snow grains across the embankment in a short-distance arrangement and give primary 
consideration to snow deposition on the windward side of the embankment in a long-distance 
arrangement. Yu et al. (2020) simulated and quantified the barrier effect of snow fence through a 
wind tunnel test and found that the higher snow fence, the better the barrier effect. The wind 
protection effect of the wind fence is improving as the fence height increases from 3.7to 7.9 m 
(Fu and Li, 2023). This study analyzes the snow-proof effect of snow fence when fence height is 
2.1, 3.0, 3.9, and 4.8 m through a single variable for the Altay—Zhundong Railway and obtains the 
best snow-proof effect when fence height is 4.8 m; this result is similar to previous results. Ma et 
al. (2022) believed that when the drifting snow flowed through the cutting, a certain scale of eddy 
currents appeared in the cutting, and the range of these eddy currents was consistent with the 
position of snow accumulation. The experimental results of Ma et al. (2022) verified the 
rationality of using a streamlined vortex to characterize the range of snow deposition in this study. 
Li et al. (2022) hold that the deposition scope of the leeward snow cover in snow fence is affected 
by fence height, which is similar to the conclusion drawn in the orthogonal tests in this study. This 
study further identified porosity as the main factor affecting the deposited snowfall on the 
pavement in the collaborative prevention and control of snow fences and embankments. In 
addition, when the arrangement distance was 20 m, the snow deposition area acted on pavement, 
and the deposited snowfall on pavement exceeded that in the absence of a snow fence, which was 
unfavorable for pavement snow protection. Therefore, when analyzing the snow cover on the 
pavement in orthogonal tests, the data obtained under the arrangement distance of 20 m were 
rejected as accidental errors. Existing embankment studies are mostly based on a single 
qualitative or quantitative analysis (Ma et al., 2021; Qiu et al., 2022). By introducing the 
mathematical statistics method, we carried out a comprehensive analysis based on single-factor 
tests and orthogonal tests in this study, and ranked the three factors in terms of their sensitivity to 
deposited snowfall in different areas. The results of this study can provide theoretical guidance 
and technical support for the prediction of snow cover distribution and the construction of snow 
protection facilities along railways. 

Railway design parameters and traffic network limit railway engineering construction. Snow 
protection measures are necessary to ensure that railway lines can successfully pass through 
extreme climate zones. This study investigates the characteristics of snow cover distribution along 
railway subgrade and the protective effect of snow fences by combining field surveys, remote 
sensing data validation, and laboratory numerical simulation. However, no laboratory wind tunnel 
tests are conducted due to the limitation of conditions. In follow-up research, wind tunnel tests 
can be used to optimize further and validate test results. Wind and snow disasters are subject to 
seasonal changes, and the erection of snow fences is limited by environment, space, cost, and 
other factors. To ensure that snow fences can serve engineering practice efficiently, stably, and 
economically, research on live snow fences may be the future trend in this field. 


5 Conclusions 


This study investigates the characteristics of snow cover distribution and the protective effect of 
different snow fences along the Altay—Zhundong Railway subgrade using single-factor tests and 
orthogonal tests. Within the study context, field survey data and global atmospheric reanalysis 
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(ERAS) climate data are utilized. The results show that in the absence of snow protection 
measures, the maximum snow cover thickness of the embankment pavement decreased by 12.6% 
relative to the cutting pavement. Hence, it is suggested that embankments be preferentially built 
in areas prone to wind and snow disasters. Setting snow fences to change the flow field 
distribution around the embankment can affect and control the number and settlement trend of 
snow particles transported on the subgrade's windward slope, road surface, and leeward slope. 
When porosity increased from 25% to 75%, the deposited snowfall on the pavement continuously 
decreased. A fence height of 4.8 m produced the largest scope of leeward snow cover and the 
smallest deposited snowfall on the pavement. An arrangement distance of 60 m produced the 
smallest deposited snowfall in front of the embankment and the smallest deposited snowfall on 
the pavement. In the collaborative prevention and control of snow fences and embankments, the 
three factors can be ranked in terms of their sensitivity to deposited snowfall within the influence 
scope of fences as follows: fence height>arrangement distance>porosity. Fence height is the main 
factor affecting the deposited snowfall within the influence scope of the off-embankment fence. 
For the embankment protection of the Altay-Zhundong Railway against wind and snow, snow 
fence with a porosity of 75%, a fence height of 4.8 m, and an arrangement distance from the 
embankment of 60 m produced the best snow control effect. 
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